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Technical  Report  Summary 


The  data  from  short-period  WWSSN  seismograms  were  used  to  study 
the  excitation  and  attenuation  of  l-Hz  crustal  phases  (P,  Pg,  Sn  and 
Lg)  in  southern  Asia.  The  data  show  that  the  attenuation  in  southern 
Asis  varies  laterally,  but  that  in  general  it  is  less  than  in  Iran  and 
the  western  United  States,  and  greater  than  in  eastern  North  America  and 
central  Asia.  The  most  severe  attenuation  occurs  in  northern  India, 
northern  Pakistan  and  Afghanistan,  where  it  is  comparable  to  the  values 
found  for  Iran  and  the  western  United  States. 

The  only  crustal  phase  with  a sharp  onset  for  southern  Asia 
earthquakes  is  P (or  Pn) . Thus  it  will  have  to  be  used  for  locating 
earthquakes  and  explosions.  For  an  = 4.0  event  this  P phase 
(vertical  component)  has  a ground  amplitude  of  at  least  1 millimicron 
to  distances  of  3500  km.  For  purely  detection  purposes,  Pg  and  Lg  will 
have  larger  amplitudes  to  distances  of  up  to  1000  km,  but  their  amplitudes 
become  less  than  those  of  P and  Sn  at  distances  beyond  1000  km  in 
southern  Asia. 

The  amplitude  variations  of  P with  epicentral  distance  in  southern 

Asia  departs  from  the  standard  Gutenberg-Richter  relation  by  as  much  as 

0.5  m,  units  at  distances  of  500  to  2500  km.  Thus  the  Gutenberg-Richter 
b 

calibration  function  should  not  be  used  at  these  distances  to  determine 
m^,  or  to  estimate  yield  of  explosions. 

From  recently  published  Soviet  studies  it  appears  that  the  attenuation 
of  seismic  crustal  phases  in  central  Asia  is  more  similar  to  that  in 


eastern  North  America  than  to  that  in  southern  Asia.  Thus  Lg  should 


be  useful  for  detecting  low-magnitude  events  out  to  distances  of  2000 
km  in  the  interior  of  the  USSR.  However,  Ur  cannot  be  used  for 
accurate  source  location,  for  which  one  will  have  to  continue  to  use 
arrival  times  of  P. 

Methods  for  inverting  higher-mode  attenuation  data  in  the  period 
range  4-10  seconds  are  being  developed  and  have  been  applied  successfully 
to  Rayleigh  wave  data  in  the  central  United  States.  The  combination 
of  fundamental-  and  higher-mode  amplitude  spectra  make  it  possible  to 
simultaneously  solve  for  source  depth  and  crustal  Q models. 

Higher-mode  data,  both  in  the  period  range  5-10  seconds  and  at 

1 second  (Lg)  have  been  combined  with  fundamental-mode  data  to  yield 

new  frequency-dependent  Q models  for  eastern  North  America.  A 

variation  of  shear  wave  Q with  frequency  (u>)  and  depth  (z)  is  assumed 

Cl 

such  that  Qq  (z,u>)  * C(z)u)  where  a is  constant.  The  results  indicate 

D 

that  the  case  a = 0 (a  frequency-independent  Q model)  do  not  explain 
the  data  very  well.  The  best  fit  to  the  entire  data  set  for  a single 
a value  for  all  periods  occurs  when  a = 0.3.  Better  fits  are  achieved 
when  a is  constant  (between  0.0  and  0.2)  over  the  period  range  5-40 
seconds  and  attain  greater  values  for  shorter  periods. 

Initial  studies  of  higher  modes  in  Iran  and  Turkey  suggest  that 
crustal  Q values  are  very  low  and  that  intrinsic  anelastlcity  is  more 
Important  than  scattering  of  seismic  surface  waves  in  the  period  range 
4-10  seconds. 
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EXCITATION  AND  ATTENUATION  OF 


SHORT-PERIOD  CRUSTAL  PHASES  IN  SOUTHERN  ASIA 


by 

Otto  W.  Nuttli 


INTRODUCTION 


Earthquakes  in  Szechwan  Province  and  Tibet,  which  occur  in  the 
border  region  separating  central  from  southern  Asia,  were  used  to 
study  the  excitation  and  attenuation  of  short-period  crustal  phases 
in  southern  Asia.  The  region  is  one  of  complex  geology  and  rugged 
topography,  where  strong  attenuation  of  short-period  crustal  phases 
and  large  variations  in  wave  amplitude  might  be  expected.  One  of 
the  purposes  of  this  study  is  to  specify  the  amplitudes  of  the  waves 
P (or  Pn),  Pg,  Sn  and  Lg  for  an  earthquake  of  a given  magnitude  as 
a function  of  epicentral  distance,  and  thus  to  determine  which  of  the 
crustal  phases  is  best  suited  for  detection  of  small-to-intermediate 
magnitude  seismic  events  in  a particular  distance  range. 

Most  of  the  data  for  the  study  came  from  the  WWSSN  seismograms 

* 

of  the  stations  KBL  (Kaboul,  Afghanistan),  QUE  (Ouetta,  Pakistan) 


and  CHG  (Chengmai,  Thailand),  for  which  we  possess  a fairly  complete 
set  of  records  for  the  years  1971-1972.  Also  used  were  data  from 
SHI  (Shiraz,  Iran),  NIL  (Nilore,  Pakistan),  KOD  (Kodaikanal,  India) 
and  SHL  (Shillong,  India),  for  which  we  have  seismograms  for  selected 
events.  Additional  seismograms  have  been  on  order  since  November 
1978,  but  they  have  not  been  delivered  as  of  April  1979.  Figure  1 
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shows  the  location  of  tlu*  earthquake  epicenters  and  of  the  stations 


used  in  the  analysts. 


METHODS  OF  DATA  REDUCTION 

All  the  earthquakes  utilized  in  the  study  were  reported  in  the 
Bulletin  of  the  International  Selsmologlcal  Centre.  The  ISC  assigned 
body-wave  magnitude  (m^>  values  to  the  earthquakes  on  the  basis  of 
teleseismic  R-wave  amplitudes,  which  are  completely  independent  of  the 
data  considered  in  this  study. 

The  amplitudes  of  T (or  Rn),  Eg,  Sn  and  Lg  were  measured  on  the 
vertical-component  seismograms,  and  of  I’g,  Sn  and  Lg  on  the 
hor Izontal-component  seismograms  as  well.  After  correction  for 
Instrument  magnification  the  amplitudes  were  equalized  to  those  of 
an  ra^  =»  5.0  earthquake  on  the  assumption  that  the  logarithm  of  the 
amplitude  of  these  waves  at  l-see  period  scales  directly  as  m^. 

Horizontal-component  amplitude  data  were  combined  with  vertical- 
component  data  on  the  basts  of  previous  observations  that  the 
resultant  horizontal  Eg  amplitude  is  approximately  equal  to  that  of 
the  vertical  component,  and  that  the  resultant  horizontal  Lg 
amplitude  is  approximately  twice  that  of  the  vertical-component. 
Inasmuch  as  the  relation  between  the  resultant  horizontal 
and  the  ver t leal-component  Sn  amplitude  is  unknown,  the  two  components 
of  the  Sn  motion  were  considered  separately. 

In  previous  studies  Nuttli  (1973,  1978)  presented  far-field  Lg 
amplitude  data  which  could  be  extrapolated  back  by  theory  to  1 km 
eplcentral  distance.  Assuming  a point-source  model  of  Lg  radiation, 
the  Lg  data  for  m^  “ 5.0  earthquakes  in  both  northeastern  America  and 
Iran  extrapolate  back  to  a value  of  800  microns  at  1 km  distance. 


L_ 


1 


J 
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This  number  should  not  bo  taken  as  a measure  of  the  Lg 'vertical- 


component  ground  motion  at  1 km  distance  for  1-sec  period  waves, 
but  rather  as  a source-related  value  which  can  be  used  to  estimate 
Lg  ground  motion  at  distances  of  50  km  and  greater.  This  can  be  done 
by  using  a theoretical  equation  for  the  attenuation  of  Lg  surface 
waves,  taking  account  of  geometrical  spreading,  dispersion  and 
anelastic  attenuation  or  absorption. 

From  the  observation  that  the  extrapolation  of  the  Lg  data  to 
l km  is  Independent  of  source  region,  it  becomes  possible  to 
estimate  the  average  value  of  the  absorption  coefficient  between  the 
epicenter  and  station  from  a single  amplitude  value.  Using  a 
family  of  curves  with  absorption  coefficient  as  parameter,  the 
average  value  of  the  absorption  coefficient  for  a path  from  the 
epicenter  to  the  station  can  be  estimated  by  plotting  the  equalized 
amplitude  (to  m^  “ 5.01  versus  epicentral  distance  and  observing  the 
curve,  corresponding  to  a particular  value  of  the  absorption 
coefficient,  on  which  the  data  point  falls. 

The  amplitude  variation  of  Pn  and  P waves  at  distances  of 
about  500  to  3500  km  does  not  follow  a simple  law,  but  rather 
depends  on  the  velocity  variation  with  depth  in  the  upper  mantle 
to  depths  of  about  1000  km.  Because  no  simple  theoretical  or 
empirical  equation  can  be  found  to  describe  the  amplitude  behavior 
at  these  distances,  t lie  amplitude  data  were  averaged  over  a 250 
km  distance  interval  and  a smooth  curve  was  drawn  through  the 
average  values  to  describe  the  variation  of  amplitude  with  distance. 
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The  amp  1 1 tudes  ot  Pg  anil  Sn  are  observed  to  decrease  mono- 
tonlcally  with  eplcentral  distance.  Kather  than  attempt  to  lit 
these  data  with  theoretical  curves.  It  Is  customary  to  tit  them  with 
empirical  straight-line  curves  on  a plot  of  the  logarithm  ot  the 
amplitude  versus  the  logarltlun  ot  the  epleentral  distance.  The 
slopes  of  the  curves  are  measures  ot  the  attenuation  ot  Pg  and  Sn, 
and  the  Intercepts  a measure  ot  the  source  excitation  tor  an 
m,  - 5.0  event. 

For  some  paths  In  southern  Asia,  1-sec  period  l.g  waves  could 
not  be  observed  on  the  seismograms,  but  waves  ot  period  1.5  to  l 
sec  could  be  seen.  In  such  cases  the  amplitudes  of  the  longer 
period  waves  were  reduced  to  those  ot  1-sec  period  waves  by  assuming 
that  the  Lg  spectrum  had  a slope  of  2 for  periods  of  l to  l sec. 

Thus  the  amplitude  of  a 1.5  sec  period  wave  was  divided  by  2.3, 
of  a 2-sec  period  wave  by  4.0  and  of  a 1-sec  period  wave  by  *) . 2 
to  reduce  the  amplitudes  to  those  of  l -sec  period  waves.  This 
procedure  can  be  justified  it  the  corner  period  for  l.g  waves  from 
earthquakes  of  m^  - 5.0  Is  greater  than  or  equal  to  1 sec.  For 
earthquakes  for  which  both  1-sec  and  longer  period  l.g  waves 
could  be  observed,  the  reduced  amplitudes  of  the  longer  period 
waves  were  found  to  be  approximately  equal  to  the  observed  amplitudes 
of  the  1-sec  period  waves,  justifying  the  assumption. 


AMPLITUDES  OK  Pn  AND  P WAVES 


Figure  2 shows  the  logarithm  of  the  vertical  component  of 
Pn  and  P waves  plotted  as  a function  of  eplcentral  distance.  Also 
shown  on  the  curve  are  the  logarthmie  average  of  the  amplitude  for 
250  km  intervals  (dashed  horizontal  lines)  and  a smooth  curve 
drawn  through  the  average  values.  The  smooth  curve  can  be  taken  to 
represent  the  average  amplitude  of  Pn  and  P (vertical  component) 
in  southern  Asia  for  an  - 5.0  earthquake.  The  scatter  in  the 
data  is  appreciable,  amounting  to  as  much  as  *0.5  m^  units  over 
most  of  the  distance  interval  of  500  to  3500  km. 

Figure  3 compares  the  average  amplitude  curve  of  Figure  2 for 
southern  Asia  with  the  Gutenberg-Richter  (1956)  amplitude  curve, 
obtained  from  their  body-wave  magnitude  calibration  function,  and 
with  amplitude-distance  data  equalized  to  m^  » 5.0  obtained  by 
Nuttli  (1971)  from  a study  of  the  amplitudes  of  P waves  from  nuclear 
explosions.  At  distances  of  500  to  2300  km  the  Gutenberg-Richter 
curve  differs  appreciably  from  the  curve  obtained  for  southern  Asia. 
At  distances  of  2300-3500  km  the  differences  are  smaller,  but  the 
Gutenberg-Richter  curve  on  the  average  tends  to  have  amplitudes  about 
0.2  to  0.3  units  larger.  The  nuclear  explosion  data,  on  the 
other  hand,  agree  remarkably  well  with  the  southern  Asia  data  at 
distances  of  1950  to  3150  km,  and  differ  by  no  more  than  0.2  m^ 
units  at  distances  of  3150  to  3500  km. 

Figure  A presents  the  P-wave  amplitude  deviations,  in  m^ 
units,  for  earthquakes  in  Tibet  and  the  surrounding  regions.  The 
number  in  parenthesis  is  the  number  of  earthquakes  used  to  obtain 
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Figure  2.  Amplitudes  of  vertical-component,  1-sec  period  P-wave  motion 
equalized,  to  an  mt>  = 5*0  earthquake.  The  dashed  lines  are  logarithmic 
averages  of  the  amplitude  data  over  250-km  intervals.  The  solid-line 
curve  is  a smoothed  curve  drawn  through  the  average  values. 
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Figure  3*  Amplitude  curves  for  vertical-component,  1-sec  period  P-wave 
motion  equalized  to  an  m-^  = 5*0  earthquake.  Iiote  the  large  departures 
between  the  Gutenberg-llichter  (1556)  curve  and  that  obtained  from  the 
data  for  southern  Asia.  The  nuclear  explosion  curve  is  from  Nuttli  (1971). 
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the  average  amplitude  deviation.  A positive  value  of  the  amplitude 
deviation  Indicates  that  the  observed  amplitudes  are  greater  than 


expected,  and  a negative  value  that  the  observed  amplitudes  are 

smaller  than  expected.  The  figure  shows  that  NIL  has  amplitudes 

2.0  times  (0.3  m.  units)  larger  than  expected,  and  KBL  lias  amplitudes 
b 

1.3  times  (0.1  m unit)  larger  than  expected.  On  the  other  hand, 
observed  amplitudes  at  SHI  and  QUF.  are  0.63  of  the  expected  value, 
and  KOD  and  CHO  0.50  of  the  expected  value.  Figure  5 presents 
similar  data  for  Szechwan  earthquakes,  for  which  the  observed 
amplitudes  at  CHG  and  QUF  are  1.3  times,  at  NIL  2.0  times  and  at 
KBL  0.79  times  the  expected  value  for  an  m^  ■ 5.0  earthquake. 

AMPLITUDES  OF  Pg  WAVES 

Figure  6 presents  the  amplitude  data  for  the  vertical  component 
of  the  Pg  ground  motion  for  southern  Asia  equalized  to  m,  - 5.0. 

The  solid-line  curve  has  a slope  of  -3.5,  Indicating  a fairly  rapid 
decrease  of  amplitude  with  distance.  By  comparison,  the  slope  of 
Pg  for  Iranian  earthquakes  recorded  In  Iran  is  -5.85  (Nutt  1 1,  1978)  at 
distances  of  700  to  1200  km.  For  the  eastern  United  States 
Evernden  (1967)  found  a slope  of  -3.42,  almost  Identical  to  the 
-3.5  found  for  southern  Asia  In  this  study.  Fvernden’s  Pg  amplitude 
values  for  the  western  United  States  are  nearly  identical  to  those 
found  for  southern  Asia.  The  Iranian  Pg  amplitudes  are  about  twice 
as  large  at  500  km  eplcentral  distance,  and  become  equal  to  the 
western  United  States  and  southern  Asia  at  1000  km.  At  larger 
distances  the  Iranian  Pg  amplitudes  decrease  much  more  rapidly. 
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A (microns) 


AMPLITUDES  OF  Sn  WAVES 


Figure  7 presents  the  amplitude  data  for  the  horizontal 

component  of  l-sec  period  Sn  waves  for  southern  Asia  equalized  to 

m,  - 5.0.  At  distances  of  750  to  900  km  the  average  curve  is 
n 

practically  Identical  to  that  given  by  Nuttli  (1978)  for  Iranian 
earthquakes  recorded  in  Iran.  However,  at  distances  beyone  900 
km  the  Iranian  curve  has  a steeper  slope,  indicating  a higher 
attenuation  in  Iran  at  the  larger  distances.  The  slope  of  the 
southern  Asia  curve  is  —3.1,  wheras  for  Iran  it  was  found  to  be 
-4.35. 

Figure  8 presents  the  amplitude  data  for  the  vertical  component 
of  l-sec  period  Sn  waves  for  southern  Asia  equalized  to  m^  “ 5.0 
Although  one  would  expect  the  average  curve  to  have  the  same  slope 
as  for  the  horizontal  Sn  data,  the  slope  of  the  vertical-component 
curve  is  smaller,  namely  -2.2.  The  differences  in  slope  result 
principally  from  differences  in  amplitude  at  distances  less  than 
1000  km,  for  which  there  are  few  data  points. 

AMPLITUDES  OF  Lg  WAVES 

Figure  9 presents  attenuation  curves  for  the  vertical  component 
of  l-sec  period  Lg  waves  for  an  - 5.0  earthquake.  The  parameter 
Is  the  absorption  coefficient,  which  is  allowed  to  vary  from  0.015  km 
to  0.0006  kjn  a range  which  spans  the  observed  values  for  short 
period  (0.1  to  1.0  sec)  Lg  waves  in  various  geological  regions.  if 
extended  back  to  l km  all  of  the  curves  would  Indicate  a ground 
amplitude  of  800  microns,  which  can  be  considered  as  the  near-field 


!■ 


\ 


! 
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A (microns) 


or  source  extrapolation  of  the  far-field  data.  In  practice,  however, 
the  curves  should  not  be  used  at  distances  less  than  about  “>0  km  to 
estimate  ground  motion,  because  1-sec  period  Lg  waves  may  not 
completely  develop  at  shorter  eplcentral  distances.  The  curves  of 
Figure  8 satisfy  the  relation 

A(A)  - 2.0855  (sin  A)-V»  (A)-1/3  exp(-JfA) 

where  A is  in  microns  and  A is  in  radians,  except  that  it  is 
customary  to  express  A in  the  exponent  in  kilometers,  and  correspondingly 
in  km-*.  The  number  2.0855  was  chosen  to  satisfy  the  far-field 
data  of  l.g-wave  amplitudes  in  both  northeastern  America  and  in  Iran, 
where  the  far-field  data  in  each  case  extrapolate  back  to  a value  of 
800  microns  at  1 km  eplcentral  distance. 

Figure  9 was  used  to  estimate  the  average  value  of  the  Lg 
absorption  coefficient  for  the  path  between  a specific  epicenter  and 
a specific  station.  The  amplitude,  equalized  to  an  m^  - 5.0  earth- 
quake, was  plotted  as  a function  of  eplcentral  distance  on  Figure  9, 
and  the  appropriate  value  of  the  coefficient  of  absorption  was  read 
from  the  curves.  The  results  of  such  analysis  are  shown  in  Figure 
10.  From  that  figure  it  can  be  seen  that  the  largest  values  for  the 
coefficient  of  absorption  occur  for  paths  between  western  Tibet  and 
KBL,  Ql'K,  and  NIL.  The  lowest  values  of  absorption  occur  for  paths 
from  eastern  Tibet  and  northeastern  India  to  NIL,  K0D  and  KBL.  By 
way  ot  comparison,  the  average  value  of  the  absorption  coefficient 
for  1-sec  Lg  waves  in  Iran  was  0.0045  km-'  (Nutt  11,  1978).  Absorption 
for  paths  across  the  Indian  shield  is  low,  but  greater  (0.0015  to 
0.0020  km-')  than  for  northeastern  America  (0.000b  km-*). 

Surprisingly,  paths  from  eastern  Tibet  and  Szechuan  across  the 
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Himalayas  to  NIL  show  relatively  small  values  for  the  absorption 
coefficient,  namely  0.0020  to  0.0025  km  * , which  Is  significantly 
lower  than  that  observed  for  Iran  and  for  the  western  United  States. 

Figure  11  presents  the  vert  leal-component  Lg  amplitudes  for 
all  the  paths  in  southern  Asia.  The  solid-line  curve  is  a theoretical 
curve  with  absorption  coefficient  equal  to  0.0025  km  If  the  curve 
is  extrapolated  back  to  1 km  it  gives  an  amplitude  of  800  microns, 
similar  to  that  found  for  Iran  and  northeastern  America.  Thus  it 
appears  that  the  source  excitation  of  1-sec  Lg  waves  is  independent 
of  source  region  for  earthquakes  with  foci  in  the  crust,  i.e.,  if 
their  m^'s  as  determined  by  teleseismic  P waves  are  equal,  their 
extrapolated  Lg  amplitudes  at  1 km  epicentral  distance  will  be 
equa  1 . 

AMPLITUDE  OF  CRUSTAL  PHASES  AS  A 

FUNCTION  OF  EPICENTRAL  DISTANCE 

For  the  determination  of  the  detection  capability  of  a seismograph 
network  it  is  necessary  to  know  which  seismic  phases  will  have  the 
largest  amplitude  at  a specific  epicentral  distance.  Figure  12 
presents  this  Information  in  graphical  form,  by  means  of  stick 
seismograms.  In  the  figure  the  amplitudes  of  the  vertical  component 
of  P,  Pg,  Sn  and  Lg  are  plotted  on  a linear  amplitude  scale  for 
distances  of  500,  1000  and  2000  km.  At  500  km  distance  in  southern 
Asia  Pg  is  the  largest  amplitude  phase,  and  Lg  is  also  comparatively 
large.  P and  Sn  are  small,  of  the  order  of  0.1  microns.  (For  an 
nc  - 4.0  event,  P and  Sn  would  have  amplitudes  of  about  10  millimicrons 


A (microns) 


A = 500  km 


A=  1000  km 


Figure  12.  Stick  seismograms  showing  the  amplitudes  of  the 

vertical-component,  1-sec  period  ground  motionfor  the  crustal 
phases  in  southern  Asia. 
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at  500  km  distance.)  At  1000  km  all  four  phases  have  about  equal 
amplitudes,  of  the  order  of  100  millimicrons  for  an  * 5.0  event. 

At  2000  km  distance  the  effects  of  absorption  on  Pg  and  Sg  are 
severe,  so  that  their  amplitudes  are  only  about  one-tenth  those  of 
P and  Sn. 

In  summary,  Pg  and  Lg  have  the  potential  for  detecting  smaller 
events  than  P and  Sn  can  detect  for  distances  up  to  about  1000  km,  but 
at  the  larger  distances  the  phases  P and  Sn  provide  the  best  detection. 
For  location  of  events,  P (or  Pn)  is  always  the  best  phase  to  use, 
because  it  is  the  only  phase  with  a sharp  onset  whose  arrival  time 
can  be  measured  to  a tenth  of  a second. 

COMPARISON  TO  AMPLITUDES  OF  CRUSTAL 
PHASES  RECORDED  IN  CENTRAL  ASIA 
Antonova  et  al  (1978)  gave  spectral  amplitude  data  for 
earthquakes  in  four  source  regions  recorded  at  the  station  Talgar, 
which  is  near  Alma  Ata  in  southeastern  Kazakh,  central  Asia.  From 
their  Figures  18,  20,  22  and  23,  which  are  attenuation  curves  for 
amplitude-to-period  ratios  of  P,  Pg,  S and  Lg  for  frequencies  of 
0.35  to  11  Hz,  equalized  to  a magnitude  5.0  event,  one  can  construct 
stick  seismograms  similar  to  those  of  Figure  12.  Figure  13  shows 
such  seismograms  for  l-Hz  waves  for  earthquake  source  regions  to  the 
northeast,  to  the  east,  to  the  south  and  to  the  west  of  Talgar.  In 
general  at  distances  of  500  and  1000  km  the  seismograms  are  similar 
in  amplitude  and  character  to  those  recorded  in  southern  Asia,  except 
that  the  Lg  waves  for  earthquakes  to  the  northeast  of  Talgar  (epicenters 
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in  Zaysan,  Altag,  Sayan,  Pribay  kalye  and  Trans-Baykal)  are  much 
larger,  of  the  character  of  what  one  would  expect  for  paths  crossing 
a shield  region,  even  though  the  region  in  question  would  be 
classified  as  tectonic  rather  than  shield.  At  2000  km  the  P waves 
for  all  four  directions  to  Talgar  are  two  to  three  times  larger  than 
for  southern  Asia,  and  again  for  earthquakes  to  the  northeast  the 
Lg  amplitudes  are  about  two  orders  of  magnitude  greater  than  those 
observed  for  southern  Asia.  Even  for  sources  to  the  east  (northeast 
China  and  Mongolia),  to  the  south  (Pakistan,  Nepal  and  India),  and 
to  the  west  (South  Tyan'  Shan',  Pamir,  Tadzhik,  Kopet-Dag,  Iran  and 
Caucasus)  the  Lg  amplitudes  at  2000  km  are  two  to  five  times  those 
observed  in  southern  Asia.  From  Figure  9 this  would  indicate  that 
the  coefficient  of  absorption  of  1-sec  period  Lg  waves  to  the 
northeast  of  Talgar  is  0.0006  km  * (similar  to  that  of  eastern  North 
America),  to  the  east  is  0.0015  km  ^ , to  the  south  is  0.002  km  * and 
to  the  west  is  0.0012  km 


CONCLUSIONS 

In  general,  crustal  phases  propagating  across  southern  Asia  are 
attenuated  more  than  in  eastern  North  America  and  in  central  Asia, 
but  less  than  in  the  western  United  States  and  in  Iran.  Only  in 
northern  India,  northern  Pakistan  and  Afghanistan  is  the  attenuation 
as  severe  as  it  is  in  the  western  United  States  and  Iran. 

P (or  Pn)  is  the  best-suited  wave  to  locate  earthquakes  in 
southern  Asia,  because  it  is  the  only  crustal  phase  with  a sharp 
onset.  It  will  have  amplitudes  of  at  least  one  millimicron  to 
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distances  as  large  as  3500  km  for  an  “ 4.0  event.  For  detection 
purposes,  Pg  and  Lg  will  have  larger  amplitudes  than  P for  distances 
to  1000  km.  Beyond  1000  km  P and  Sn  have  the  largest  amplitudes. 

The  P amplitudes  for  southern  Asian  earthquakes  depart  by  as 
much  as  0.5  m^  units  or  greater  from  the  Gutenberg-Richter  (1956) 
amplitudes  at  distances  of  500  to  2500  km.  Tlius  at  these  distances 
the  Gutenberg-Richter  calibration  function  should  not  be  vised  to 
determine  m^  of  southern  Asia  earthquakes.  The  nuclear  explosion  data 
given  by  Nuttll  (1971)  provide  a close  fit  to  the  southern  Asia  data, 
but  unfortunately  the  explosion  data  curve  only  begins  at  2200  km 
distance . 

From  data  given  by  Antonova  e t a 1 (1978)  the  attenuation  of 
1-sec  period  Lg  waves  for  earthquakes  in  Baykal  recorded  at  Talgar 
(43°N,  78°E)  is  similar  to  that  in  eastern  North  America,  even  though 
the  former  is  for  a path  across  a tectonic  region  (ust  north  of 
Mongolia.  If  a broad  generalization  can  be  made  from  this  one  data 
set,  one  might  expect  that  similar  small  values  of  l.g  attenuation 
would  be  found  in  the  Russian  platform  and  the  non-tectonic  Siberian 
Interior.  Thus  Lg  could  well  have  the  potential  for  detecting  events 
of  magnitude  even  less  than  4 out  to  distances  of  2000  km  in  the 
interior  of  the  USSR,  although  Lg  could  not  be  used  to  accurately 


locate  such  events. 


REFERENCES 


Antonova,  L.V.,  F.K.  Aptieaev,  R.J.  Curochkina,  I.L.  Nersesov 
A.V,  Nikolaev,  A, I,  Ruzaykln,  E.N.  Sedova,  A.V.  Sitnikov, 
l.S.  Tregub,  L.D.  Fedorscaya  and  V.I,  Halturin  ( 1 *» 78) . 
Experimental  Science  Studies  of  the  Earth's  Interior, 

Institute  of  Earth  Physics,  USSR  Academy  of  Sciences,  Nauka, 
Moscow.  * 


Evernden,  J.l.  (19p7).  "Magnitude  determination  at 
near-regional  distances  in  the  United  States, 
Soc.  Am..  57,  591-619. 


regional  and 
Bull.  Sels. 


Gutenberg,  B.  and  C.F.  Richter  (1956).  "Magnitude  and  energy  of 
earthquakes,"  Ann.  Geo f Is.  (Rome).  9,  1-15. 

Nuttll,  O.W.  (1971).  "The  amplitudes  of  teleseismic  P waves  " 
Bull.  Seism.  Soc.  Am..  61,  343-356. 


Nuttll,  O.W.  (1973).  "Seismic  wave  attenuation  and  magnitude  relations 
ot  eastern  North  America,"  Jour.  Geoph.  Res..  78,  876-885. 

Nuttll  O.W.  (1978).  "Amplitudes  and  attenuation  of  seismic  crustal 
phases  at  regional  distances  in  Iran"  in  Research  in  Seismolovv. 

' R<;Porl.  bV  B- J-  Mitchell  and  O.W.  Nuttll,  ARPA  Contract  ^ 
144620-73-C-0042,  Nov.  1978,  12-31. 


Determination  of  Q models  from  Surface  Wave  Spectra 

by 

Chiung-chuan  Cheng 

The  determination  of  Q structure  from  surface  wave  spectra  of 
a single  seismogram  (our  previous  reports)  requires  the  source 
mechanism  as  well  as  the  focal  depth  be  known.  These  are  not  always 
known  as  precisely  as  we  would  like.  However,  Tsai  and  Aki  (1970) 
Investigated  variations  of  the  surface  wave  amplitude  spectrum  which  were 
due  to  uncertainties  of  these  factors  and  concluded  that  effects 
caused  by  factors  other  than  the  focal  depth  were  relatively  small. 

We  have  made  further  comparisons  between  the  nature  of  spectral 
variations  caused  by  varying  focal  depth  and  that  by  varying  Q 
structures. 

A stochastic  Inversion  scheme  is  applied  to  obtain  Q from 
surface  wave  spectra. 
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The  effect  of  focal  depth  and  Q structure  on 
surface  wave  spectra 

The  Rayleigh  or  Love  wave  spectral  amplitude  of  a certain 
mode  at  a certain  period  is  assumed  to  be 

V - y0e~YR  (1) 

where  y0  is  computed  amplitude  without  attenuation. 

The  effect  of  focal  depth  is  present  in  yQ  through  a multipli- 
cation factor,  the  excitation  coefficient,  which  in  turn  is  the  source 
spectral  amplitude  multiplied  by  a very  complicated  function  of  the 
focal  depth  and  the  source  mechanism. 

The  wave  amplitude  decreases  exponentially  with  the  product  of 
the  attenuation  coefficient  y and  the  epicentral  distance  R.  For  a 
layered  earth  model,  the  Rayleigh  wave  attenuation  coefficient  is 
related  to  the  compress lonal-  and  shear-wave  velocities,  a and  tf, 
and  the  associated  quality  factors,  and  (Anderson  et_  aj_. , 1965; 
Mitchell,  1975). 

"N  S£_\  _!  N ( ^L_  \ 

Y " Tti'cr  3<V  wp3  Qcti+^C2  3MwpaVJ  (2) 

The  complexity  of  the  way  surface  wave  spectrum  varies  with  the 

focal  depth  is  illustrated  in  figures  1 and  2.  The  contours  show  the 

values  of  log  y.  The  source  is  assumed  to  be  a 45°  thrust  fault  having 

24 

a seismic  moment  of  1.0  x 10  dyne-cm.  The  epicentral  distance  is 

1000  km,  and  the  station  is  located  at  an  azimuth  of  45°  from  the  source. 

A Q model  which  has  Q0  values  of  250  from  the  surface  to  a depth  of 

n 

17  km,  2000  at  all  greater  depths,  and  a being  twice  everywhere 
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(Herrmann  and  Mitchell,  1975)  Is  used.  Figure  1 shows  the  result  of 
the  fundamental  mode,  and  figure  2 the  first  higher  mode  Rayleigh 
wave.  The  general  feature  shows  that  the  energy  lobes  migrate  to 
the  longer  periods  as  the  focal  depth  increases.  Nevertheless,  their 
relative  heights  vary. 

Figure  3 shows  the  effect  of  varying  the  Q values  in  the  upper 
layer.  The  same  source  is  assumed  and  is  fixed  at  a depth  of  10 
km.  The  contours  show  log  y for  the  fundamental  mode  Rayleigh  wave. 
From  equation  (1),  it  is  seen  that  log  y is  a function  linear  in  the 
varying  Q values.  A straight  line  drawn  for  constant  period  in  the 
figure  will  intersect  the  contours  with  equ i- intervals,  showing  a 
constant  slope  for  log  y to  decrease  with  increasing  • 


Stochastic  inversion  of  surface  wave  spectra 


An  observed  spectral  amplitude  z is  in  general  different  from 
the  theoretical  amplitude  y computed  from  equation  (1),  with  presumed 
seismic  moment  and  Q model.  To  fit  y and  z,  the  seismic  moment  and 
Q model  have  to  be  adjusted. 

Assuming  that  values  are  twice  that  of  everywhere,  from 
equations  (1)  and  (2), 


Rtt 

T 


N 

L 

f.=  l 


[■»( 


a£_) 

3V 


utpB  + 


A q 


-l 


+ RC 


C.n  y - 2.n  z 


where  AQ^ J is  the  adjustment  for  in  each  layer  of  the  model,  and 
the  seismic  moment  is  adjusted  to  by  multiplied  by  exp  (-RC) . 


Fitting  y and  z over  all  the  periods  and  modes  of  available 
data  points  yields  a system  of  equations,  which  can  bo  written  as 
Ax  - b 

where  x is  the  vector  to  be  solved, 

x " (AQsr  a^b2*  A(^$3*  •••  A(Vn’  c^T* 

With  the  transformations  (Wiggins,  1972)  that 


A' 

- s“s 

: 

A W 

x’ 

- W~S 

X 

b* 

- s"*5 

b 

tile  matrix  equation  can  be  rewritten  as 
A'x'  - b' 

W and  S are  assumed  to  be  diagonal  matrices.  The  stochastic  inversion 
is  applied  to  solve  for  x'. 

x'  - (A’ V + 0 1)  A ' p ' 

Thus, 

v,  T - 1 2-1  -V,  t - 1 

x - (W  A S A + o i ) W ^ A S b (3) 

This  Inversion  scheme  yields  an  estimate  of  x that  minimizes 

e2  * iA'x'  - b')T  ( A ’ x ’ - b’)  + a2  x’V 

- <S-S  Ax  - S"S  b)T  (S'*5  Ax  - S-^  b)  + 

O2  (W  ‘Sc)1'  (W  ‘Sc) 

The  introduction  of  S enables  us  to  add  weights  to  the  observations. 
I lie  combination  of  o and  W provides  a measure  of  how  freely  each  of  the 
parameters  can  deviate  from  its  presumed  value. 

Spectral  comparison  for  vertical  component  Rayleigh  waves  from 
the  southeastern  Missouri  earthquake  of  October  21,  1965  (Mitchell,  1973; 


' 
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Herrmann,  1974)  recorded  at  MDS  was  made  with  a two-layer  Q model  (our 
previous  report).  A discrepancy  could  be  noticed  for  the  fundamental 
mode  spectrum  at  short  periods  (figure  4),  so  the  observed  data  set  has 
been  Inverted  using  equation  (3)  In  a hope  to  Improve  the  fit.  The 
starting  Q model  Is  the  same  as  the  one  used  in  spectral  comparisons 
before.  However,  the  upper  layer,  which  was  1 7 km  thick  and  had  Q 

n 

values  of  250,  Is  divided  Into  a layer  of  7 km  overlying  another  of 
10  km;  and  a layer  of  50  km  Is  added  to  the  top  of  the  half-space, 
which  had  a of  2000. 

The  resulting  model  Is  shown  In  table  1.  A comparison  of 
observed  and  theoretical  spectra  appears  In  figure  5. 

Discussion 

The  way  In  which  the  surface  wave  spectrum  varies  with  the  focal 
depth  is  different  from  the  way  It  varies  with  the  Q model,  and  is 
much  more  complicated.  Attempts  to  Invert  surface  wave  spectra  to 
obtain  the  focal  depth  and  Q model  simultaneously  is  impractical 
compared  to  that  to  obtain  the  Q model  only,  because  of  Its  non- 
linearity and  large  computer  time  consumpt ion.  Even  if  the  focal 
depth  Is  not  well  determined  in  advance,  the  Q model  can  still  be 
obtained  by  inverting  with  several  presumed  focal  depths.  An  adequate 
fit  between  the  theoretical  and  observational  spectra  occurs  with  the 
proper  focal  depth.  Analyses  of  several  sets  of  spectral  data  using 
this  method  are  in  progress. 
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Figure  2.  Contours  showing  the  first  higher  mode 
Rayleigh  wave  logarithmic  spectrum  for 
various  focal  depths. 
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PERIOD  (SEC) 


Figure  3 


Contours  for  the  fundamental  mode  Rayleigh 
wave  logarithmic  spectrum  attenuated  with 
a two-layer  model  which  has  various  values 
in  the  upper  layer. 


AMP  (CM-SEC) 


MDS  Z 21  OCT  65 


Figure  4.  Comparison  of  theoretical  and  observational 
amplitude  spectra,  based  on  a two-layer  Q 
model. 
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Source  Studies  of  Some  Eurasian  Earthquakes 

by 

Hui-Vuin  Wen  Wang 


In  the  last  semi-annual  report,  we  mentioned  that  in  order  to 
know  the  source  time  functions  of  the  earthquakes,  several  determinations 
should  be  made  for  comparing  synthetic  seismograms  with  observed  ones: 

(1)  source  velocity  model;  (2)  receiver  velocity  model;  (3)  seismic 
wave  attenuation  data;  (4)  fault  plane  solution;  (5)  focal  depth; 

(6)  seismic  moment.  Especially  the  seismic  velocity  model  and 
attenuation  coefficients  are  the  two  most  important  effects  on 
synthetic  seismograms.  In  this  study,  surface  waves  generated  by  22 
earthquakes  which  occurred  in  the  Tibetan  Plateau  and/or  on  the  Tibetan 
border  were  used  to  analyze  the  crustal  structure  across  Tibet. 

Geologic,  seslmic,  and  gravity  data  imply  that  the  current 
tectonic  activity  in  Asia  is  the  consequence  of  continental  collision 
between  India  and  Eurasia.  East-west  trending  strike-slip  faults  are 
most  observed  in  China.  Movements  on  them  may  allow  material  lying 
between  the  stable  portions  of  the  Indian  and  Eurasian  plates  to  move 
laterally  out  of  the  way  of  those  two  plates.  The  Tibetan  Plateau  has 
a unique  topographic  feature  of  high  altitude,  which  is  not  a direct 
result  of  the  continental  collision  but  of  uplift  during  northerly 
under thrust ing  of  the  Indian  subcontinent  along  the  deep  crustal 
fracture.  This  led  to  an  approximate  doubling  of  crustal  thickness 
under  the  Tibetan  Plateau  (Molnar  & Tapponmier,  1975;  Powell  & Conaghan, 
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1973).  Santo  (1965),  Gupta  6 Naraub  (1967),  Chon  6 Molitar  (1975),  and 
Chun  & Yoshli  (1977)  analyzed  the  dispersion  of  the  fundamental  Rayleigh 
waves  and  found  that  the  crustal  thickness  In  Tibet  was  65  to  70  km. 

Table  1 gives  the  events  which  are  used  In  the  study.  The  paths 
of  use  in  the  construction  of  the  Tibetan  velocity  model  are  obtained  by 
selecting  those  which  are  largely  confined  to  the  Tibetan  region. 
Therefore,  the  surface  waves  to  stations  as  KHl.,  QllE,  Nil.,  ND1,  and  SHL 
have  been  digitized  and  proceed  by  the  multiple-filter  method  (Oziewonski 
£t  al . , 1969).  Thirty  paths  are  analyzed  to  obtain  the  group  velocities 
of  fundamental-mode  Rayleigh  waves  in  the  6-60  second  period  range. 

The  averaged  values  and  standard  deviations  are  given  in  Table  2.  The 
crustal  structure  was  perturbed  by  trial  and  error  until  the  theoretical 
dispersion  curve  best  matched  the  observed  data.  First,  we  tried  Chun 
& Yoshii  (1977)'s  TP-4  model.  Their  curve  (triangular  sign  in  Figure  1) 
is  a quite  good  fit  to  our  data  (solid  line)  for  periods  smaller  than 
30  seconds,  whereas  the  group  velocities  are  relatively  lower  at  longer 
periods.  Also  TP-4  possesses  a trough  which  is  somewhat  broader  and 
shifted  in  period  when  compared  with  our  data,  which  have  a minimum 
value  centered  at  30  seconds.  Chun  & Yoshii's  model  is  characterized 
by  a crust  of  about  70  km  in  thickness  with  low  shear  velocities  and 
the  upper  mantle  velocity  is  low  with  shear  velocity  of  about  4.45  km/sec. 
The  parameters  of  the  TP-4  model  are  given  in  Table  3.  By  increasing 
the  shear  velocity  of  the  lower  crust  to  3.81  kin/see  (named  model  A in 
Table  3),  this  new  model  indicated  by  crosses  in  Figure  1 shows  the  best 
fit  to  the  observed  data  at  longer  periods.  It  tends  to  possess  slightly 
lower  velocity  values  for  periods  smaller  than  30  seconds,  the  largest 
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difference  being  0.02  km/sec.  Such  differences  are  negligible  and 
cannot  be  resolved  within  the  accuracy  limits  of  the  observed  data. 
Because  the  seismic  stations  we  chose  are  near  to  the  Tibetan 


region,  the  structure  complexities  may  not  affect  the  surface  waves 
greatly.  It  is  not  unreasonable  to  use  the  same  velocity  model  for 
the  receiver  crust  and  upper  mantle  structure. 

Determination  of  attenuation  coefficients  for  paths  across  Eurasia 
is  very  difficult  because  the  paths  between  sources  and  stations  usually 
cross  more  than  one  tectonic  province.  Yacoub  & Mitchell  (1977) 
divided  Eurasia  into  two  regions,  one  considered  to  be  relatively  stable 
and  the  other  considered  to  be  tectonic  in  nature.  By  employing 
Mitchell's  method  (1975)  to  fit  observed  to  theoretical  amplitudes  by 
least  squares  for  determining  Rayleigh-wave  fundamental-mode  attenuation 
coefficients,  Yacoub  & Mitchell  (1977)  tested  four  Eurasian  earthquakes 
and  found  that  the  tectonic  regions  exhibit  higher  attenuation  than 
the  stable  regions  in  the  period  range  below  about  20  sec,  whereas  in 
the  period  range  above  20  sec,  no  clear  difference  can  be  observed  for 
the  two  regions.  We  can  take  their  tectonic  attenuation  data  and 
apply  inversion  theory  (Backus  & Gilbert,  1970)  to  obtain  models  of 

Q * as  a function  of  depth  beneath  Tibetan  region.  After  Q*  models 

p p 

are  available,  theoretical  surface  wave  attenuation  coefficient  values 
for  all  modes  can  be  computed  by  using  the  formula  of  Anderson  £t  a.i 
(1965). 

The  collection  of  70  mm  film  chips  of  WWSSN  stations  at  the 
department  is  not  complete  enough  to  allow  us  to  make  fault  plane 
solution  and  focal  depth  determination  for  any  one  of  the  Tibetan 
earthquakes  listed  in  Table  1.  Therefore  we  used  ISC  data  and  plotted 


the  first  P-wave  motions  on  a lower  focal  sphere.  The* most  likely 
solutions  are  shown  in  Figure  2.  They  are  all  equivalent  to  a 
combination  of  reverse  and  strike-slip  faulting,  with  onenodel  plane 
dipping  east-west  and  the  other  north-south.  The  T-axis  trends  southwest. 

Discussion 

Molner  et_  a_l . (1973)  showed  that  the  strike  directions  of  most 
earthquakes  in  Tibet  trend  either  north-south  or  east-west  with  a 
northeast  trend  of  the  P-axis.  They  are  normal  faulting  with  an 
eastj-west  T-axis  whichmay  reflect  east-west  flow  of  material  in  the 
lower  crust  and  upper  mantle  beneath  Tibet  to  compensate  for  the 
pressure  imposed  by  the  plate  motion.  The  apparent  contradiction  we 
found,  that  the  Tibetan  is  under  compression,  may  be  resolved  by 
observing  more  data  by  ourselves.  We  are  planning  to  make  a trip  to 
the  National  Geophysical  Data  Center  in  Boulder  in  May. 

As  Figure  1 shows,  although  the  velocity  model  we  computed 
has  a relatively  higher  group  velocity  at  longer  periods  as  compared 
with  TP-4,  the  latter  data  are  still  within  our  confidence  limits. 


I 
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Figure  1,  Dispersion  curves  for  Tibetan  Plateau.  Crosses 

represent  velocities  for  model  A.  Triangles  represent 
group  velocities  for  model  TP-4.  Solid  line  indicates 
the  observed  data  from  30  paths  across  Tibetan 
Plateau  with  standard  deviations  marked  by  vertical 
lines. 
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Table  3 


Tibetan  velocity  model 


model 

thickness 

(km) 

P velocity 
(km/sec) 

S velocity 
(km/ sec) 

densll 

(R/cm 

TP-4 

3.50 

4.50 

2.60 

2.40 

8.50 

5.98 

3.45 

2.80 

16.00 

5.98 

3.42 

2.80 

10.00 

5.80 

3.37 

2.75 

30.00 

6.30 

3.64 

2.90 

7.70 

4.45 

3.30 

A 

3.50 

4.50 

2.60 

2.40 

8.50 

5.98 

3.45 

2.80 

16.00 

5.98 

3.42 

2.80 

10.00 

5.80 

3.37 

2.75 

30.00 

6.60 

3.81 

2.90 

7.70 

4.45 

3.30 
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Frequency  Dependence  of  in  the 
Continental  Crust 


Brian  J.  Mitchell 


Introduction 


Inversions  of  surface  wave  attenuation  data  to  obtain  internal 
friction  (Q^S  models  of  the  crust  and  upper  mantle  usually  assume  that 
the  degree  of  frequency  dependence  of  Q0  is  small  enough  to  be  ignored. 

p 

Most  models  of  the  continental  crust  and  upper  mantle  (e.g.  Mitchell, 

p 

1973b;  Lee  and  Solomon,  1975;  Herrmann  and  Mitchell,  1975;  Mitchell,  1975) 
are  consequently  dependent  upon  the  assumption  that  Q0  is  independent 

p 

of  frequency.  The  purpose  of  this  report  will  be  to  combine  higher  and 
fundamental  modes  to  test  the  appropriateness  of  the  frequency-independence 
assumption  in  the  inversion  of  seismic  surface  wave  attenuation  data. 

The  phase  Lg  is  often  the  largest  phase  on  short-period  seismograms 
recorded  at  regional  distances.  In  the  frequency  band  10  Hz  to  about 
2 seconds,  Lg  consists  of  a superposition  of  numerous  higher-mode 
Rayleigh  waves  (Knopoff  ejt  al^. , 1973).  At  longer  periods  (about  4- 
10  s) , it  is  often  possible  to  isolate  one  or  two  higher  modes.  If 
the  nature  of  the  frequency  dependence  of  Qa  can  be  established,  then 

p 

it  should  be  possible  to  predict  the  attenuation  of  high-frequency 
(1  Hz)  Lg  waves  from  longer-period  higher-mode  data,  and  vice  versa. 


Data 


P 


The  most  extensive  set  of  surface  wave  attenuation  data  at 
intermediate  periods  in  continental  regions  has  been  obtained  in  eastern 
North  America  (Herrmann  and  Mitchell,  1975).  The  only  higher  mode 
attenuation  data  in  the  period  range  4-10  seconds  have  also  been 
obtained  for  that  region  (Mitchell,  1973a;  Herrmann,  1973).  Fundamental- 
mode  attenuation  data  in  the  period  range  4-40  s will  be  combined  with 
data  at  a period  of  0.4  s (Kurz  and  Redpath,  1968)  to  obtain  new  Qg^ 
models  of  the  crust  in  eastern  North  America.  The  models  which  are 
obtained  are  then  used  to  calculate  higher-mode  attenuation  coefficients 
for  comparison  with  the  observed  attenuation  coefficients  of  the  higher 
modes  between  4 and  10  s and  with  the  observed  attenuation  coefficients 
for  Lg  at  shorter  periods  in  eastern  North  America  (Nuttli,  1973; 

Street,  1976;  Bollinger,  1979;  Pomeroy,  1977;  Nuttli  and  Dwyer,  1978). 

The  details  of  the  method  are  described  below. 


Procedure  and  Results 

The  amplitude  of  particle  motion  of  the  higher  Rayleigh  modes, 
even  at  relatively  short  periods,  is  significant  to  large  depths  in 
the  crust  (Figure  1).  Amplitudes  of  the  higher  modes,  for  this  reason, 
can  be  compared  to  those  of  the  fundamental  mode  to  investigate  the 
frequency  dependence  of  shear  wave  internal  friction  (QgS  in  the  crust. 

The  first  step  in  the  process  of  investigating  the  frequency 
dependence  of  Qg  in  the  crust  is  to  obtain  models  of  Qg*  as  a function 
of  depth  from  fundamental-mode  Rayleigh-wave  attenuation  data.  In 


T 


H 
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Normalized  Amplitude 


obtaining  those  models  I assume  that  Qg  is  frequency  dependent  such 
that  Qg(o>, z)  - C(z)u>a  where  w is  angular  frequency  and  a is  a constant. 
The  case  a - 0 corresponds  to  the  case  where  Qg  is  independent  of 
frequency.  It  is  possible  to  satisfy  the  fundamental  mode  data  with  a 
wide  range  of  values  for  a.  Two  cases  (a  = 0.0  and  0.5)  appear  in 
Figure  2. 

The  next  step  is  to  do  the  forward  problem  of  calculating 
attenuation  coefficients  for  the  higher  modes  using  the  models  obtained 
from  inversion  of  the  fundamental-mode  data.  A difficulty  in  this 
process  is  that  the  observed  higher-mode  data  sometimes  correspond  to 
a single  mode  and  sometimes  to  several  modes.  This  problem  is 
illustrated  by  Figures  3 and  4 where  theoretical  group  velocities  for 
the  eastern  United  States  model  of  Mitchell  and  Herrmann  (1979)  and 
theoretical  spectral  amplitudes  for  that  model  and  selected  event- 
station  examples  are  plotted.  It  is  apparent  that  higher-mode 
observations  at  velocities  of  about  3.5  km/sec  in  the  eastern  United 
States  consist  predominantly  of  the  first  higher  mode  at  periods 
greater  than  about  3 seconds.  At  shorter  periods  there  is  much  overlap 
in  group  velocities  for  the  various  modes  and  several  modes  contribute 
to  the  observed  Lg  wave  form.  However,  Lg  attenuation  coefficients  are 
usually  determined  for  waves  travelling  at  group  velocities  between 
3.3  and  3.6  km/sec.  If  that  is  the  case,  an  average  attenuation 
coefficient  for  the  first  four  higher  modes  should  approximate  the 
attenuation  coefficient  for  Lg  at  a period  of  1 second.  At  a period  of 
2 seconds,  the  first  two  higher  modes  should  be  sufficient. 

Figure  5 presents  observed  higher-mode  attenuation  coefficients 
and  theoretical  coefficients  for  various  values  of  a,  assuming  that 
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Observed  attenuation  coefficients  for  the  fundamental  Rayleigh  mode. 
Open  circles  and  vertical  lines  are  the  values  and  standard  deviations 
observed  by  Herrmann  and  Mitchell  (1975).  The  closed  circle  is  a data 
point  from  Kurz  and  Redpath  (1968). 


Attenuation  Coefficient  (X IO"3),  km*1 


a is  a constant  throughout  the  entire  period  range  between  1 and  40 
seconds.  It  is  clear  that  a frequency  independent  (a  = 0.0)  Qg*  model 
will  not  satisfy  the  higher-mode  attenuation  data.  Further,  it  appears 
that  none  of  the  models  for  which  a is  a constant  throughout  the  entire 
period  range  fit  the  data  very  well.  The  case  a = 0.3,  perhaps,  comes 
closest  to  explaining  all  of  the  data. 

Additional  models  were  explored  in  which  a was  allowed  to  vary 
with  frequency.  Attenuation  coefficients  for  some  of  these  models 
appear  in  Figures  6 and  7.  A very  good  fit  to  the  data  was  achieved  for 
the  case  in  which  a was  set  equal  to  zero  at  periods  greater  than  4 
seconds  and  equal  to  0.5  at  periods  of  4 seconds  or  less.  The  hump  in 
the  curve  corresponds  to  the  jump  in  a at  a period  of  4 seconds  and  can 
be  smoothed  out  with  a less  abrupt  change  in  a.  Good  fits  to  the  data 
can  also  be  achieved  by  letting  a increase  gradually  with  decreasing 
period  at  periods  less  than  4 seconds.  Two  such  cases  appear  in  Figure 
7 where  a is  set  at  values  of  0.0  and  0.2  at  long  periods.  Undoubtedly, 
several  other  distributions  of  a over  the  period  range  1-40  seconds  will 
also  satisfy  the  higher-mode  attenuation  data. 

There  are  several  Qg*  models  for  the  crust  of  eastern  North  America 
which  might  explain  the  fundamental  and  higher-mode  data.  A frequency- 
independent  model  and  a model  for  the  case  of  constant  a (=  0.3)  over 
the  entire  frequency  range  between  1 and  40  seconds  appear  in  Figure  8. 
Some  models  for  which  a varies  with  period  are  presented  in  Figure  9. 
Comparisons  of  the  theoretical  attenuation  coefficient  values  to  the 
observations  (Figures  5-7)  indicate  that  the  cases  for  which  Ot  increases 
with  decreasing  period  are  better  than  those  for  which  a is  a constant. 
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coot  I lclents  to  observed  values  for  various 
values  of  i*.  (1)  a - 0.0  for  T N 4 seconds, 

o - 0. S for  1 < T < 4 seconds.  (2)  a - 0.0 
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seconds.  (1)  n - 0.0  for  T > 4 seconds,  0.4 
v a ' 0,0  for  1 v T v 4 seconds. 
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Higher-Mode  Attenuation  in  the  Middle  East 

Iran  to  Turkey 
by 

Brian  J.  Mitchell 

Introduction 

Methods  previously  developed  to  determine  shear  wave  internal 
friction  models  in  North  America  are  being  used  to  study  the  attenuation 
of  higher-mode  Rayleigh  waves  in  the  Middle  East.  Although  this  work  is 
still  In  progress,  preliminary  results  have  been  obtained  which  are 
important  to  understanding  surface  wave  propagation  and  attenuation 
through  a portion  of  the  Middle  East  from  eastern  Iran  to  western 
Turkey. 


Data 

As  described  in  other  sections  of  this  report,  higher  modes  with 
periods  of  a few  seconds  will  be  combined  with  fundamental-mode  data 
to  study  the  anelastic  properties  of  the  earth's  crust.  The  method 
requires  the  use  of  sources  for  which  the  depths  and  fault-plane 
solutions  are  known.  Fault-plane  solutions  for  several  suitably 
located  earthquakes  in  the  Middle  East  have  been  determined  by  McKenzie 
(1972).  They  are  listed  in  Table  1 and  appear  on  a map  in  Figure  1. 

The  paths  cover  an  east-west  trending  region  from  eastern  Iran  to 
western  Turkey,  most  of  which  is  comprised  of  plateaus  and  small 
mountain  ranges. 

The  essential  data  consist  of  Rayleigh  wave  amplitude  data, 
both  fundamental  and  higher  modes.  In  addition  to  the  amplitude  data. 
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however,  group  velocity  data  are  also  required  in  order  to  establish 
a crustal  model  for  that  region.  Both  these  types  of  data  have  been 
obtained  using  the  multiple-filter  method.  The  group  velocities  are 
plotted  in  Figure  2.  It  is  interesting  to  note  that  the  higher-mode 
velocities  are  observed  only  to  periods  as  great  as  about  6 seconds.  By 
comparison,  higher  modes  are  observed  in  eastern  North  America  to 
periods  as  large  as  about  10  seconds.  The  scatter  in  both  the 
fundamental-  and  higher-mode  Rayleigh  wave  group  velocity  data  is  quite 
large,  probably  because  of  the  complexity  of  the  crust  in  that  region. 

The  group  velocities  are  quite  low,  indicating  that  the  shear  velocities 
in  the  crust  of  that  region  of  the  Middle  East  are  lower  than  those 
which  would  be  expected  for  a shield  or  stable  region. 

Method  and  Preliminary  Results 

In  order  to  apply  our  method  to  obtain  attenuation  data  and  Q 
models  for  the  portion  of  the  Middle  East  between  eastern  Iran  and 
western  Turkey,  it  is  first  necessary  to  obtain  a model  of  seismic 
velocities  and  densities.  Since,  to  our  knowledge,  no  crustal  model 
is  available  for  that  region,  it  is  necessary  to  obtain  a model  by 
using  the  group  velocity  data  of  the  present  study.  Because  of  the 
relatively  large  scatter  in  the  observations,  it  was  decided  to 
obtain  a model  by  simple  tr ial-and-error  fitting  of  the  data.  The 
solid  lines  in  Figure  2 indicate  the  theoretical  group  velocities  for 
the  model  in  Table  2.  It  is  interesting  to  note  that  the  model  is 
very  similar  to  the  model  of  Bucher  and  Smith  (1971)  for  the  Colorado 
Plateau  in  the  western  United  States,  differing  only  by  having  one 
layer  which  is  l*j  km  thicker  than  the  corresponding  layer  in  the  Colorado 
Plateau. 
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Figure  2.  Observed  fundamental-mode  (closed  circles)  and  higher-mode  (open  circles)  Rayleigh  wave 
group  velocities  across  Iran  and  Turkey.  The  solid  line  indicates  theoretical  values 
for  the  crustal  model  in  Table  2. 
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The  availability  of  a crustal  model  and  sources  wfth  known 
depths  and  fault-plane  solutions  permits  the  calculation  of  surface 
wave  amplitude  spectra  (Levshin  and  Yarson,  1971;  programmed  by 
R.B.  Herrmann)  for  the  fundamental  and  higher  modes.  Examples  of 
these  spectra  appear  in  Figures  3 and  4,  where  they  are  compared  with 
observed  spectra  for  the  earthquake  of  27  July  1966  assuming  two 
different  focal  depths.  Figure  3 compares  observed  spectral  amplitudes 
with  theoretical  values  corresponding  to  a source  depth  of  40  km  and 
Figure  4 compares  the  observed  values  to  theoretical  values  for  a 
source  depth  of  20  km.  The  greater  depth  is  the  value  given  by  the 

I USGS  location.  In  both  cases  the  observed  higher  mode  spectral 

amplitudes  are  much  smaller  than  those  predicted  by  theory.  This 
result  indicates  severe  attenuation  of  higher  modes,  either  by  low  Q 
values  in  the  crust  or  by  scattering  produced  by  inhomogeneities.  If 
the  focal  depth  for  the  earthquake  is  40  km,  it  is  likely  that  scattering 
of  the  higher  modes  is  most  important,  since  the  fundamental  mode 
spectral  amplitudes  in  Figure  3 exhibit  little  attenuation.  If  the 
focal  depth  is  near  20  km,  however,  anelastic  attenuation  is  likely 
to  be  the  important  mechanism,  since  both  higher  modes  and  the  short-period 
fundamental  mode  spectral  amplitudes  in  Figure  4 exhibit  large 
attenuation.  The  position  of  the  spectral  minimum  in  the  observations 
at  a period  of  about  40  seconds  corresponds  more  closely  to  the  theoretical 
spectral  hole  for  the  case  of  the  shallower  source.  This  preliminary 
result  indicates  that  low  Q is  likely  to  characterize  the  crust  of  Iran 
and  Turkey.  Further  studies  should  provide  enough  data  to  yield  a Q 
model  for  that  region  which  will  explain  both  fundamental-mode  .and 


H 


\ 


i 


M 


I j 


-78- 


higher-mode  attenuation  data.  Frequency-dependent  models  will  be 
sought  which  will  also  explain  the  attenuation  of  regional  phases 

such  as  Lg. 
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